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Ab initio results for structural and electronic properties 
of NaCl-type FeN are presented. Calculations are performed 
within a plane- wave and ultrasoft pseudopotential framework. 
Competition among paramagnetic, ferromagnetic and some 
possible antiferromagnetic phases is examined. We find the 
ferromagnetic phase stable overall as observed. Stabilization 
over the unpolarized phase is obtained by splitting one flat 
feg-type band crossing the Fermi energy. Comparison with 
CrN reveals that the magnetic ordering of CrN, which consists 
of double- ferromagnetic sheets compensating along the [110] 
direction of the cubic cell and stabilized by the distortion of 
the squared symmetry on (100) plane, is not effective for FeN. 
The different behavior of FeN and CrN can be traced back to 
the smaller ionicity and magnetization of the former. 
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I. INTRODUCTION 

Stable compounds^ formed by Ee and N, such as 7'- 
Fe4N and C-Fe2N, or quasi-stableB like a"-FeigN have 
been known for some time. On the other hand, al- 
most nothing-jwas reported about FeN until the work of 
Heiman et alB, where for the first time, they synthesized 
Fe x Ni_a; with 0.4 < x < 0.75. In this work FeN was 
found to be ferromagnetic. Successive studies found FeN 
in NaCl-type crystaL-structurea, and non-ferromagnetic 
at room temperatures. More recent works.are in contra- 
diction with the earlier ones: Morita et alB reported FeN 
antifcrromagnet at temperatures below 100 K. Moreover, 
Suzuki et al.U determined the FeN structure to be of ZnS 
fee type. They also found their sample to exhibit a mic- 
tomagnetic character that was suggested be related to 
the antiferromagnetism of FeN. 

Regarding the structure, it seems clear that FeN can 
assume both zincblende and rocksalt structure, depend- 
ing on the sample history, and both of them may contain 
some deficiencies of N atoms. Regarding magnetic prop- 
erties, recent Mossbauer spectroscopy measurements sug- 
gested that NaCl-type FeN is antiferromagnetictjEI, with 
Neel temperature around room temperature. However, 
Mossbauer spectroscopy is not able to determine com- 
pletely magnetic ordering. 

For both structural and magnetic properties ab initio 
calculations are suited to furnish a prediction about the 
ground state of FeN, and to help guide understanding 
of experimental results. Here we investigate the rela- 
tive stability of paramagnetic (PM), ferromagnetic (FM) 
and some possible antiferromagnetic (AFM) phases of 
NaCl-type FeN that have been observed or suggested 



for other transition-metal nitride (or oxide or carbide) 
compounds. Compared to other transition-metal nitrides 
(TMNs), it has been noticedB that NaCl-type FeN shows 
electronic properties close to those of CrN that we re- 
cently investigatedlj, e.g. the Fermi energy falls exactly 
on a pronounced peak of the density of states, mainly 
formed by Fe d-like flat bands. In Ref. [ll] the authors 
come to the conclusion, on the basis of the rigid band 
model to evaluate the DOS and the Stoner condition of 
ferromagnetism, that only CrN, FeN and CoN in rocksalt 
phase are likely to be stabilized by magnetic ordering. 

Comparison with CrN is particularly interesting be- 
cause it undergoes a cubic-to-orthorhombic transition ac- 
companying the PM to AFM ordering at a Neel temper- 
ature close to the room temperature, with AFM phase 
formed by double ferromagnetic sheets compensating 
along the [110] direction. The structural distortion is 
decisive in determining the most stable among the dif- 
ferent AFM phases, and stress relief is likely to be the 
ultimate driving force towards distortion. If such an in- 
volved mechanism of magnetic-ordering is driven by the 
high DOS at Fermi level, it could be argued that a similar 
magnetic ordering occurs for FeN as well. Thus, for pos- 
sible AFM phases, we decide to consider here the [110]- 
double and single sheet AFM arrangements (indicated 
as AFM^ 110 j a nd AFM^ 110 j, respectively and described 

fully in Ref.lMj), as well as the AFM[m] phase (very com- 
mon for NaCl-type transition metal oxides), consisting 
of single ferromagnetic sheets alternating along the [111] 
direction. 

Equilibrium lattice constants, energies, and magnetic 
moments are presented for all the considered phases. 
Also, band structures and density of states (DOS's) are 
presented for PM and FM phases. Our local-spin-density 
calculations are performed iH a plane- wave and ultrasoft 
pseudopotential frameworklia Use of ultrasoft pseudopo- 
tentials allow us to obtain well converged results for a 
cutoff energy equal to 30 Ryd. We used sets of 10 to 110 
special k-points (depending on the structure) for the self- 
consistent calculations, up to 280 k-points to evaluate the 
DOS. 



II. STRUCTURE 



Experimentally it was found for the lattice constapt of 
FeN in rocksalt and zincblende phase Oq s =8.5 bohra and 
ag b =8.1-8.2 bohrEH3, respectively. This result is some- 
what surprising; indeed, assuming the same anion-cation 
bond length in the two phases, the relation between lat- 
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though this crude approximation considerably overesti- 
mates the actual difference between the two structures, 
it is generally true that in most cases the rocksalt phase 
have larger volume. Doubts about the measurements are 
confirmed by a recent theoretical workliil where struc- 
tural properties for transition metal nitrides are calcu- 
lated. They found for a value in good agreement 
with experiments, but a much smaller a™ ~ 7.5 bohr in 
the paramagnetic phase. The discrepancy with experi- 
ment is so large that it seems unlikely it could be due to 
different magnetic orderings. Also, in Rcf. |ll| the equi- 
librium lattice constant for other phases (ferromagnetic 
and a couple of antiferromagnetic phases) was evaluated, 
but no relevant changes on ag s were found. Finally, the 



value i 



7.5 bohr is consistent with that of other rock- 
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salt TMN structures. For instance, experiments assign to 
AFM CrN a value a™ ~ 7.8 bohr, and we should expect 
for FeN a somewhat smaller lattice constant. 

To try to explain the disagreement with experiments 
we consider here other AFM phases. In particular 
among the AFM phases examined in Ref. 
AFM[ 110 ] was not taken into account. The 
AFM^ 110 j transition is found to produce, for CrN, a large 
increase of lattice constant!^ (~ 2%), although not large 
enough to account for the discrepancy with experiments 
of FeN, and the largest energy gain with respect to the 
PM phase. In Table | we report our calculated values of 



the 
PM to 



~ 2 meV). 

Thus, notwithstanding the similarity of DOS at Fermi 
level, CrN and FeN are stabilized by very different mag- 
netic order and structure. In Table [| we report the mag- 
netization m (per formula unit, i.e. per couple Fe-N) 
for the considered phases of FeN. The FM phase has the 
highest magnetization (1.65 /is), of which ~ 0.15 (j,b 
comes from polarization of N, and ~ 1.5 fis from polar- 
ization of Fe. For CrN the magnetization is more that 2 
fiB per formula unit. Also, the charge transfer is differ- 
ent. For FeN it is equal to about one electron, so that 
the approximate ionic configuration Fc +1 N _1 can be de- 
duced, i.e. Fe is near a d 7 configuration. In case of CrN, 
the charge transfer is larger (close to 2 electrons), due 
to the higher electronegativity of Fe with respect to Cr. 
i.e. to the smaller electronegativity difference with N. 
The higher covalency of FeN contributes to it having a 
smaller lattice constant, whereas CrN has larger ionicity 
and polarization. 

In Rcf. [l(] it is argued that for CrN the energy gain 
from the planar distortion is due to a relief of tensile 
stress stored in the bonds between metal atoms with 
antiparallel spins. This stress can be conveniently de- 
fined as the excess stress produced by a transition from 
the AFM^ 110 j phase (where there is no stress by symme- 
try) and the AFM^ 110 j cubic phase with the same lattice 
constant. The transition produces a planar anisotropic 



for the corresponding magnetic orderings. Evidently, stress, whose component T fl y 01 is tensile, i.e. the spin 



contrary to what happens for CrN, the structure of FeN 
is basically unaffected by magnetic ordering, and only a 

0. 4% increase of (Iq s is found in the FM phase with re- 
spect to the PM. Our results are in substantial agreement 
with Ref. [ll] (they found a =7.47 and 7.54 for PM and 
AFM phases, respectively). 

From our calculations, it results that NaCl-typc FeN 
is ferromagnetic. The energy gain of FM phase with re- 
spect to the PM is less than 0.1 eV per formula unit, 

1. e. significantly smaller than in case of CrN. Also, the 
competition between FM and the examined AFM phases 
is very close (energy differences are few hundredths of 
eV). To be confident about the precision of our calcula- 
tions, we evaluated the FM-AFM energy differences by 
performing FM calculations in the same symmetry of all 
the AFM structures considered, and tested the differ- 
ences for sets of special K-points of increasing size. We 
conclude that, among the AFM phases, the lowest in en- 
ergy is AFM^ 110 j , also the lowest among the cubic phases 

of CrN0. In that case, the orthorhombic [110] shear 
distortion of CrN produces a further energy gain for the 
AFM^ 110 j structure (anisotropic in the (100) plane due to 
magnetic ordering) that causes the latter to be the most 
stable phase overall. To investigate the effect of planar 
distortion on FeN we applied the same distortion exper- 
imentally observed for CrN (i.e. a ~2% reduction of the 
bisection angle at the squared base of the tetragonal cell. 
Surprisingly, we find very little change in the energy with 
respect to the cubic AFM^ 110 j phase (the increase is just 



[no] 

antiparallel metal atoms tend to get closer each other. In 
light of our results for FeN, we speculate that this mag- 
netic order-derived stress is sensitive to the magnitude 
of magnetization. Of course, the different bond length 
(slightly larger for CrN than for FeN) plays a role as well 
and could justify in itself a different behavior. 



III. ELECTRONIC PROPERTIES 

In Figure the DOS for the PM phase of FeN is shown. 
In this energy region, the total DOS (top panel) is almost 
completely made up by d states of Fe and p states of N 
(bottom panel). The d-like DOS shows two high peaks 
close to the Fermi energy (Ep). The lowest in energy 
is completely occupied and is composed by a mixture of 
75% t2 g -type and 25% e g -type character-. It is this peak 
which falls at in paramagnetic CrNEJ (a rigid band 
model works well for the PM phase of this system). The 
second peak in the d DOS is entirely t2 g in character, 
and is centered at Ef. The peaks located between -8 
eV and -4 eV are mostly N p, with a mixture of ti g - 
type and (predominantly) e ff -type states. The N p-like 
DOS gives an appreciable contribution to the DOS at the 
Fermi energy as well, unlike CrN for which only d-type 
states contribute. This DOS peak at Ep arises from a 
strong dpa antibonding band. 

The effect of FM ordering is easily visible in Figure 
^. Now the Fe d-like DOS (middle panel) is shown de- 
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composed in £ 2g and e g states. The two d-type peaks 
are now split by an exchange splitting of roughly 1.6 eV. 
After this splitting both peaks in the spin up channels 
are completely occupied, whereas the upper peak of the 
spin down channel is empty. Thus, 3 electrons per Fe 
atom saturate almost completely three t\ g states, and ~ 
1.5 electron charge of mostly ti g type remains unpaired. 
Consistently with a cation to anion charge transfer of ~ 
1 electron, i.e. with a d 7 configuration for Fe, there are 
2.5 electrons occuping nearly paired e g states. The N p 
character in the DOS contributes to the magnetization, 
reflecting the N moment of ~ 0.15 \ib- 

The same features can be seen in the band structure. 
In Figure || the PM FeN band structure is presented. It 
is easy to distinguish the very flat i2g-type band crossing 
the Fermi level through A, Z and Q directions, that pro- 
duces the half-occupied peak of PM DOS shown in Figure 
|]. After FM ordering (Figure [I]), the band is split, in such 
a way that the up band (left panel) is shifted ~ 1.5-2 eV 
down to the Fermi energy, and the down band rises up 
slightly, but still remains partially occupied. Evidently, 
much of the minority DOS peak at ~0.8 eV arises from 
regions away from the symmetry lines plotted in Figure 



A feature observable in Figure [| is that the spin split- 
ting for some bands is strongly k-dependent, unlike the 
Stoner picture in which the splitting is rigid. For ex- 
ample, all bands at T are clustered within 1 eV in the 
paramagnetic case and remain so in the minority bands, 
whereas in the majority bands they span a range of nearly 
2 eV. Such non-rigid behavior has been observed in other 
transition metal compounds, and reflects a combination 
of (i) an anisotropic exchange potential (around the Fe 
atom) and (ii) hybridization between d orbitals on a 
strongly magnetic atom (Fe) and p orbitals on a weakly 
magnetic atom (N). 



IV. CONCLUSIONS 

In this paper the structure of NaCl-type FeN has been 
investigated by ab initio calculations. The ferromagnetic 
phase is found to be stable against the paramagnetic and 
all the antiferromagnetic phases here considered. How- 
ever, the energy differences between the magnetic phases 
are within few hundredths of eV, and this fact, coupled 
with sensitivity of the magnetism to stoichiometry, could 
explain the experimental variations reported in the lit- 
erature. In order to get insights into the magnetic or- 
dering mechanisms, a comparison with CrN was made. 
Although for both compounds a sharp o?-like peak is 
present at Ej?, the differences in charge transfer, magneti- 
zation and differing wavefunction character at Ep lead to 
very different magnetic properties. In particular, the or- 
thorhombic distortion that is able to stabilize AFM CrN 
has no effect on FeN. Since the distortion indicates the 
presence of tensile stress stored on bonds between metal 



atoms with antiparallel spins, we suspect that the stress 
could depend sensitively on the magnitude of magnetiza- 
tion. 

The relief of stress is now accepted to be an important 
factor in the structure and properties of surfaces and 
interfaces. Our results strongly suggest that the stress 
relief may be instrumental in determining the magnetic 
order in compounds such as FeN. 
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TABLE I. Equilibrium lattice constant do (in bohr), en- 
ergy E (in eV/formula unit) and magnetic moment m (in 
/is/metal atom) of NaCl-type FeN calculated for different 
magnetic phases. Energies are referred to the lowest one, i.e. 
that of FM phase. 
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FIG. 1. Top panel: total and atom by atom decomposed 
DOS of PM NaCl-type FeN. Bottom panel: Fe d-type and N 
p-type contributions to the total DOS. In the shown energy 
window this two angular components are largely dominant. 
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FIG. 2. Top panel: total DOS of FM NaCl-type FeN for 
spin up and down charge. Middle: Fe d £2 9 -type (solid 
line) and e s -type (dashed) contributions to DOS. Bottom: 
N p-type contribution. 
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FIG. 3. Band structure for PM NaCl-type FeN. A flat band 
of t2 3 character running from F to L almost in touch with the 
Fermi level is responsible for the pronounced peak of DOS at 
the Fermi Energy. 
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FIG. 4. Band structure for FM NaCl-type FeN. Left and 
right panels refer to spin up and down bands, respectively. 
Directions and K-points are the same of the previous Figure. 



•5 



